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ABSTRACT 



We use arcminute-scale data from the Sunyaev-Zel'dovich Infrared 
^sO | Experiment to set limits on anisotropies in the cosmic microwave background 

radiation in open and spatially-flat-A cold dark matter cosmogonies. There are 
no 2-cr detections for the models tested. The upper limits obtained are consistent 
with the amplitude of anisotropy detected by the COBE/DMR experiment. 
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1. Introduction 

With the detection of anisotropy in the cosmic microwave background (CMB) by the 
COBE/DMR experiment (Smoot et al. 1992; Wright et al. 1992; Bennett et al. 1996) 
and a host of other detections on angular scales greater than ~ 10' (Ganga et al. 1994; 
Hancock et al. 1996; Piccirillo et al. 1997; Netterfield et al. 1997; Gundersen et al. 1995; 
Tucker et al. 1997; Piatt et al. 1996; Masi et al. 1996; Tanaka et al. 1996; Inman et al. 1996; 
Griffin et al. 1997; Scott et al. 1996), CMB anisotropy measurements have become a 
powerful tool for testing models of cosmic structure formation. At smaller angular scales, 
however, the detection of anisotropies has proven more difficult (Griffin et al. 1997; 
Myers et al. 1993; Subrahmanyan et al. 1993). Moreover, calculation of theoretical CMB 
anisotropies at these angular scales is difficult and has only recently been done accurately. 
Thus, the detection of arcminute scale anisotropies and their classification is one of the 
significant hurdles to be overcome in the characterization of the CMB. 

In this paper we use arcminute-scale CMB anisotropy data from the Sunyaev-Zel'dovich 
Infrared Experiment (SuZIE), an experiment designed to measure the Sunyaev-Zel'dovich 
(SZ) effect in distant clusters of galaxies. These data were taken as part of a program 
to characterize the instrument for SZ observations, but can be used to set interesting 
limits on the power spectrum of CMB anisotropies. Theoretical predictions for primary 
CMB anisotropy spectra are used in conjunction with these data to set limits on the 
normalization of cosmic structure formation models in low-density open and spatially-flat 
(with a cosmological constant A) cold dark matter (CDM) universes. 



2. Data 

SuZIE was originally conceived and constructed to make millimeter- wave observations of 
both the thermal and kinetic SZ effects in clusters of galaxies from the Caltech Submillimeter 
Observatory on Mauna Kea (Holzapfel et al. 1997a; Wilbanks et al. 1994). It employs 
a 2x3 array of bolometers cooled to 300 mK, with beam widths of approximately 1.7' 
(FWHM) and beam throws of ±2.3' in right ascension (see insets of Figure [IJ also Figure 2 
of Church et al. 1997). The instrument is described fully by Holzapfel et al. (1997b). 

Each pair of detectors in a single row is electronically differenced in an AC bridge, 
with the output of each bridge phase-synchronously demodulated to produce a DC signal 
proportional to the difference in power absorbed by the two detectors (Rieke et al. 1989; 
Wilbanks et al. 1990; Glezer, Lange, & Wilbanks 1992). In terms of atmospheric rejection, 
this is equivalent to a square wave chop between the two detectors at infinite frequency. 
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Observations were made at 140 GHz (~ 11% bandpass) and were performed by fixing the 
telescope and using the rotation of the Earth to make the sky 'drift' (in right ascension) 
through the array's view. This eliminates spurious signals arising from mechanical chopping 
or motion of the telescope. 

The output from the AC bridge, after calibration, gives the difference in temperature 
between two spots on the sky, 

A tj =Tt- Tj, (1) 

where T, and Tj are the CMB temperature at the two spots on the sky. The focal plane 
has two rows of three pixels. The rows are separated by 2.0' in declination, with the 
three pixels in each row separated by 2.3' in right ascension. The data used here are the 
differences between the two outermost pixels in each row, detectors 3 and 1 in row 1 and 
detectors 6 and 4 in row 2, which are separated by 4.6' (Figure 1; also see Figure 2 of 
Church et al. 1997). 

Two 36' regions of the sky were observed, one centered at R.A. 10 h 21 m 49.00 s and Dec. 
4° 4' 23.00" (Region 1) and the other at R.A. 16 h 30 m 17.00 s and Dec. 5° 56' 0.00" (Region 
2; epoch 1950). The instrument is Nyquist sampled at 5 Hz. The data have an offset and 
linear drift removed, and are coadded and binned into 45" bins. The data are shown in 
Figure [I] and the zero- lag window function associated with one row of pixels, Wi, is shown 
in Figure |2|. The uncertainty in the gaussian beam widths, estimated at 5%, translate to 
errors in limits on models primarily through the uncertainty introduced into the calibration. 
The beam width uncertainties are ignored otherwise (Ganga et al. 1997). The calibration, 
obtained from scans of Uranus, is known to 8%. The data reduction and calibration are 
described more fully in Church et al. (1997). 



3. Models 

The CMB fractional temperature perturbation, 5T/T, can be decomposed into a sum 
over spherical harmonics according to 

Srp OO I 

~r(M) = £ E a lm Y lm (0A)- (2) 

l=2m=-l 

In this scheme, a convenient characterization of the CMB anisotropy is its angular power 
spectrum, Ci, defined in terms of the ensemble average 

\ a lmO'ii m ') — Ci5w5 mm i. (3) 
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Prior to the construction of semi-realistic cosmogonies, it was conventional to use a 
gaussian parameterization for the CMB anisotropy spectrum. The fiducial CDM (fCDM) 
model, an Einstein-de Sitter cosmological model with gaussian, adiabatic, scale-invariant 
energy-density perturbations (Harrison 1970; Peebles & Yu 1970; Zel'dovich 1972), focussed 
attention on what has now come to be known as the flat angular spectrum (Peebles 1982), 

6C 2 24vr (Q/Tp) 2 

' 1(1 + 1) 5 1(1 + 1)' { ' 

In this context, Q, sometimes referred to as Q r ms-PS; is the implied quadrupole amplitude 
of the CMB anisotropy spectrum and To is the present CMB temperature. Q should not be 
confused with the actual amplitude of the quadrupole moment of the CMB anisotropy, to 
which SuZIE is not sensitive. This spectrum is shown as the straight, solid line in Figure ^. 

The flat spectrum is simply a large-scale approximation of the fCDM spectrum. On 
smaller scales, the pressure of the photon-baryon fluid causes the spectrum to deviate from 
this approximation. The fCDM spectrum, with standard recombination, h = 0.5, and 
baryonic-mass density parameter Qb = 0.0125/i -2 , is shown as the dotted line in Figure 0. 
Here, the Hubble parameter H = 100/t km s _1 Mpc -1 . Although still of some historical 
interest, the fCDM model, which has a clustered-mass density parameter Q = 1, is now 
known to be an inadequate representation of the observed universe. 

Low-density CDM cosmogonies, with flat spatial sections and a cosmological constant 
or with open spatial sections and no A, are consistent with a large fraction of present 
observations. For fiat- A models, see Stompor, Gorski, & Banday (1995), Ostriker & 
Steinhardt (1995), Ratra et al. (1997), Liddle et al. (1996), and Ganga, Ratra, & 
Sugiyama (1996), and for open models see Kamionkowski et al. (1994), Ratra et al. 
(1997), Ganga et al. (1996), and Gorski et al. (1996). Since the low-density A-CDM 
models are assumed to have flat spatial sections, the simplest power spectrum for gaussian, 
adiabatic, energy- density perturbations is the scale-invariant one. The simplest spectrum 
for these perturbations consistent with open spatial sections is that which is generated 
by quantum-mechanical fluctuations during an early epoch of inflation in an open model 
(Ratra & Peebles 1994, 1995; Bucher, Goldhaber, & Turok 1995; Yamamoto, Sasaki, & 
Tanaka 1995). The values of the parameters Qq, h and Qb which characterize the spectra 
in the low-density CDM models used here are chosen to be roughly consistent with present 
observational estimates of f2 , h, the age of the universe, and the constraints on Qb that 
follow from the observed light element abundances in the standard nucleosynthesis model 
(Ratra et al. 1997; see Table 1 for parameter values). A standard recombination open model 
spectrum and a standard recombination flat-A model spectrum are shown in Figure |2|. 

The general procedure for the computation of CMB anisotropy spectra is discussed 
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by Sugiyama (1995). Since SuZIE is sensitive to sub-arcminute-scale anisotropies, we have 
computed the open model spectra to I = 7500, where photon diffusion strongly suppresses 
the power. Photon diffusion cuts off power on larger scales in the fiducial and flat- A CDM 
models, so we have only computed the spectra in these models up to I — 3500. 



4. Analysis 

As a starting point for our analysis, we note that Church et al. (1997) have shown 
that the reduced x 2 ' s f° r each of the data sets used here are less than one. This gives 
us a good indication that the data does not show a detection; i.e., it is consistent with 
experimental noise. As shown below, this is indeed the case. To set numerical limits on 
various cosmological models, however, and in order to compare the results obtained for 
different models, we use a more refined technique. 

A primary goal of CMB anisotropy observations is to measure the angular spectrum 
of the fluctuations. Given the width of the SuZIE window function, it is not possible to 
extract a value for each multipole moment to which the experiment is sensitive. Rather, one 
can hope, at best, to determine a few observational numbers with reasonable accuracy. It is 
therefore necessary to simplify the problem by assuming a model for the spectrum over the 
range of I to which the experiment is sensitive. Typically, one assumes a functional form 
for the shape of the spectrum as a function of /, with the overall normalization allowed to 
be a free parameter. The shape of the spectrum can also depend on other free parameters. 
One then compares this family of spectra to the data and determines the value of the 
normalization and other parameters which best reproduce the data. 

For the purpose of the following discussion we assume that the SuZIE sky signal is 
purely CMB without any offset or gradient, which would have been removed during data 
reduction. Following Bond et al. (1991) to account for the baseline removal, the likelihood, 
at the nominal calibration, corresponding to a given spectrum is 

L oc . X e~* 2/2 , (5) 

v /det(/ T M- 1 /)det(M) 

where, 

X 2 = A t (m _1 - M- 1 f(f T M- 1 f)- 1 f T M~ lT )A. (6) 

Here, A is the vector of temperature differences and M is the correlation matrix of the 
data in the context of the model considered. The correlation matrix, M, is the sum of 
correlations arising from experimental noise, S, and those from correlations intrinsic to 
the CMB, C. That is, M = C + S. In the case of SuZIE, the noise correlation matrix 
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is approximately diagonal, with the diagonal elements simply equal to the variance of the 
corresponding measurement. This is true so long as atmospheric correlations and other 
correlations induced in the process of reducing the data can be ignored (Church et al. 
1997). C is calculated according to 

tf« = i£(2I + l)CiWtf, (7) 



47T , 



where 



167T 






sin(m$ fe /2) n '" 
m$ fe /2 



[1 - cos (m$ )] Pim (Oi) Pim (Oj) cos [m {(pi - 4>j)} 



Here Bu = exp[— 1(1 + l)of /2] is the Legendre transform of the i th beam, assumed to be 
gaussian, $^ is the extent of a single data bin in right ascension, 0.75' x cos 5 in this case, 
and $o is the azimuthal separation of the two points differenced in a single sample, 4.6'. 
6 is the declination of the pixel above or below the center of the array (±1'), and P/ m is 
defined by 

Y lm (0A) = Pim(0)-e im t (9) 

See Bond (1996) and White & Srednicki (1995) for general discussions of CMB window 
functions. In eqs. @ and @, / is the array of functions that have been fit out of the data 
(in this case an offset and gradient from each channel; it is thus an 96 x 4 matrix for each 
region) . 

With this prescription, all models yield the same likelihood at Q = /xK (i.e., nothing 
but noise), providing a convenient way to normalize across models. 

Following Ganga et al. (1997), we account for calibration uncertainty by taking 

£ = -=J x f°° dQ' e-W^ 2 ^^L(Q'), (10) 

V27T a c Q "' 

where ac is the fractional uncertainty in the calibration, 0.08. This is conservative since it 
assumes the calibration uncertainty is completely correlated in all detectors. This is the 
case if there is an error in the model of the calibration source, but does not apply if the 
uncertainties are statistical, due to, for example, uncertainty in the beam measurements. 

This calibration-uncertainty-corrected likelihood is computed for all the models 
considered. Assuming a uniform prior in Q gives a posterior probability density distribution 
equal to the likelihood function. 

For the models considered here there are no 2-<r anisotropy detections, defined 
according to the prescription discussed in Ganga et al. (1997; also see Berger 1985). To 
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quantify the constraints, we quote two different 2-cr upper limits for each model and data 
set considered. The highest posterior density (HPD) prescription provides the smallest 
possible 2-0" upper limit. In this scheme, the upper and lower 2-cr limits on Q, Q u and Qi, 
are defined such that 



£ ■ AQ = 0.9545 x / C-dQ, (11) 

Jo 

and Q u — Q\ is minimized. Here, Qi is zero for all models, indicating that there are no 
detections. The equal tail (ET) prescription gives a larger 2-cr upper limit, calculated using 

£ -dQ = 0.9772 x / C-dQ. (12) 

o Jo 

Finally, we note that expected sky rms can be written as a sum of contributions from 
different multipoles, ST rms 2 = J2i(^T rms 2 )i, where 

{STmJ^ = T'^-^-QWl (13) 

47T 

In Figure [| we plot this quantity for the four representative spectra shown in Figure |2]. This 
plot indicates that though the SuZIE window function peaks at I ~ 2400, with l c g ~ 2340 
and Z e -o.5 ~ 1330 and 3670 (as defined by Bond 1996), the SuZIE sensitivity is model 
dependent. In these models SuZIE is most sensitive to multipoles I ^ 1000, and is still 
sensitive down to I ~ 200. 



5. Results and Discussion 

Table 1 lists the 2-cr upper limits on Q. Also listed are the limits on the 
bandtemperature, 5Ti, obtained using the upper limits on Q, equation (|T3"D, and 



ST, 



\£ 



1=2 



"""' (14) 



1(1+1) 



Note that 5TJ is not a function of I (Bond 1996, also see Ganga et al. 1997, eq. [7]). 
Like <5T rms , 5Ti should be fairly independent of the model, but unlike ST rms , it can also 
be compared to results from experiments with different window functions. The degree to 
which 5Ti varies from model to model gives an indication of how good the flat bandpower 
approximation is. 

Figure f| shows representative likelihood functions for the combined SuZIE data and 
some of the models considered. 



In all cases, the SuZIE 2-er upper limits on Q are consistent with, and less restrictive 
than, the detections obtained by the COBE/DMR experiment. In general, these limits are 
lower for the open models than for the fiat-A models. The 2-<r upper limits on STi are much 
less model dependent than the Q limits. There remains, however, a ~ 25% difference in STi 
between the two extreme cases, indicating that using the flat bandpower approximation to 
calculate the expected signal that should be seen by an experiment is not accurate at these 
angular scales. 

Like open CDM cosmogonies, anti-tilted flat-A models (n > 1) predict more power 
on arcminute angular scales than do flat-A models without tilt (n = 1). The SuZIE 2-cr 
upper limits on Q for models with Q = 0.3, h = 0.6, QbIi 2 = 0.0175 and n = 1.2 and 1.4 
are 18 [iK and 11 [iK (HPD), respectively. These can be compared to the limits for the 
twenty-first model of Table 1, which has n — 1. When these results are combined with the 
COBE/DMR result, it is likely that spectral indices of order 1.5 or larger can be ruled out 
for this particular cosmogony. 

Atmospheric noise is a concern for all ground based CMB anisotropy experiments. 
Because we have assumed that the experimental noise is gaussian, if it were to be 
significantly non-gaussian it could change the limits set here. In order to check this, we 
have examined histograms of the data in each bin and found them sufficiently gaussian. We 
have also repeated the above analysis including estimates of the correlations induced by the 
atmosphere (Church et al. 1997). In all cases, these limits are comparable to or less than 
those quoted above without the atmospheric corrections. Contributions from dust at these 
angular scales are not expected to be significant (Tegmark & Efstathiou 1996). 

Since the telescope is generally stationary while data is taken, SuZIE is mostly immune 
to problems caused by motion of the instrument. Between scans, however, the telescope 
slews, raising the possibility of temperature drifts in the instrument at the beginning of a 
new scan. To check that this is not an issue, we repeated the above analysis excluding the 
first seven bins in each scan. The results were qualitatively the same as those from the 
entire data set, with limits a bit less strict, numerically consistent with the loss of data. 
Church et al. (1997) have done a more sophisticated analysis by correlating the data with 
the temperature of the SuZIE 300 mK stage and come to the same conclusion. 

We have checked that the truncated spectra used here does not affect the results by 
"filling in" the C\ spectra to / = 10000 using the C\ value at the largest computed I value. 
As one might expect from an inspection of Figure |3], the numbers do not change to the 
accuracy quoted here. 
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6. Conclusion 

The SuZIE has been used to set limits on primary anisotropies in the CMB at 
arcminute angular scales. Though there are no detections in this data set, the likelihood 
functions shown in Figure [| peak tantalizingly close to the maximum likelihood Q values 
obtained from the COBE/DMR data, lending hope that CMB fluctuations will soon be 
detected on arcminute scales. Even if they are not, however, increased sensitivity may allow 
for meaningful limits to be set on cosmological parameters such as Qq or the spectral index 
of primordial fluctuations. 

The next generation version of the SuZIE instrument, SuZIE 1.5, is now operational 
and has taken data simultaneously in three wavelength bands: 1.1, 1.4, and 2.1 mm 
(Mauskopf et al. 1997). These three bands correspond to three atmospheric windows near 
the peak brightness of the CMB and close to the minimum in the confusion expected from 
astrophysical sources on these angular scales (Fischer & Lange 1993). The new receiver's 
ability to observe simultaneously in these three bands will allow residual atmospheric noise, 
which limited the sensitivity of the previous instrument, to be identified by its frequency 
correlation between channels and then removed. 

We thank Thor Wilbanks for contributions to the SuZIE instrument and for help with 
the observations. This work was supported in part by NASA grant NAGW-4623 and by 
National Science Foundation grant AST-95-03226. BR was supported in part by National 
Science Foundation grant EPS-9550487 and matching support from the state of Kansas. 
NS was supported by Grant-in-aid for Scientific Research No. 8740193 from the Ministry of 
Education of Japan. The CSO is operated by the California Institute of Technology under 
funding from the National Science Foundation, contract AST-93-13929. 
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Table 1: SuZIE Region 1, 2, and Combined 2-cr Upper Limits on Q and STi (in /xK) for 
Various Models a 
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a The first 13 models are open CDM, the 14 th is fCDM, the 11 following fCDM are flat- A CDM models and, 
finally, the flat spectrum model is at the bottom. The quoted values of STi are determined using the data 
from both Regions 1 & 2. 

fc The HPD prescription is described in eq. [ O . 
c The EQ prescription is described in eq. |l2|. 
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Figure Captions 

Fig. 1. — Measured temperature differences on the sky as a function of scan position. The 
upper two panels show data from the two detector rows for Region 1, the lower two panels 
show data from the two detector rows for Region 2. As stated in the text, each of the 
four data sets correspond to the differences between two of the six pixels, as represented 
schematically by the six circles in each plot. Each pixel is separated from it's neighbors 
by 2.3' in right ascension and 2' in declination. The filled circles indicate the pixels being 
differenced for each data plot. See Church et al. (1997) for a full description of the data 
reduction. 

Fig. 2. — SuZIE zero-lag window function (Wj = Wm in eq. |J; scale on right axis) and 
four primary CMB anisotropy spectra normalized to the DMR two-year data (scale on left 
axis). The open, flat-A and fCDM spectra are computed for Qsh 2 = 0.0125. The open 
CDM, fCDM and flat spectra are normalized to the DMR galactic-frame maps (Gorski et 
al. 1994, 1995; also see Ratra et al. 1997), the flat-A spectrum is normalized to the DMR 
ecliptic- frame maps (Bunn & Sugiyama 1995; also see Ratra et al. 1997), and in all cases 
the quadrupole moment is excluded from the analysis of the DMR data. Since the DMR 
four-year data have been completely analyzed for only some of the models considered here, 
we use the DMR two-year data normalizations, bearing in mind that the four-year data 
indicates values for Q which are ~ 10% smaller than those obtained with the two-year data 
(e.g., Gorski et al. 1995, 1996). The flat bandpower l-o (HPD) upper limit obtained here is 
placed at the location of the window function peak. See, however, Figure |3[ 

Fig. 3. — (ST rms 2 )i (eq. [T^]), as a function of I for the CMB anisotropy spectra shown in 
Figure |2|. These curves should be compared to the SuZIE window function shown in Figure ^. 
Note the multiple "sensitivity" peaks for each spectrum, and that the major contributions to 
the signals arise from multipole values lower than those indicated by considering the window 
function (Figure 2) alone. 

Fig. 4. — Likelihood functions given the combined Regions 1 and 2 SuZIE data, for the four 
CMB anisotropy angular spectra of Figure 0. 
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